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Cosmic bubble collisions provide an important possible observational window on the dynamics of 
eternal inflation. In eternal inflation, our observable universe is contained in one of many bubbles 
formed from an inflating metastable vacuum. The collision between bubbles can leave a detectable 
imprint on the cosmic microwave background radiation. Although phenomenological models of the 
observational signature have been proposed, to make the theory fully predictive one must determine 
the bubble collision spacetime, and thus the cosmological observables, from a scalar field theory 
giving rise to eternal infiation. Because of the intrinsically non-linear nature of the bubbles and 
their collision, this requires a numerical treatment incorporating General Relativity. In this paper, 
we present results from numerical simulations of bubble collisions in full General Relativity. These 
simulations allow us to accurately determine the outcome of bubble collisions, and examine their 
effect on the cosmology inside a bubble universe. We confirm the validity of a number of approxi- 
mations used in previous analytic work, and identify qualitatively new features of bubble collision 
$H spacetimes. Both vacuum bubbles and bubbles containing a realistic inflationary cosmology are 

studied. We identify the constraints on the scalar fleld potential that must be satisfied in order to 
obtain collisions that are consistent with our observed cosmology, yet leave detectable signatures. 



^ I. INTRODUCTION 

in 

O Accelerated expansion of the universe and spontaneous symmetry breaking, two powerful ideas in modern physics, 

^ conspire in many theories to give rise to a phenomenon known as eternal inflation. This phenomenon occurs when 
a region of the universe is in a metastable vacuum with positive energy density. Classically, these vacua are stable, 
fS| and drive inflation (exponential expansion). Quantum mechanically, they can decay via the nucleation of expanding 
J> bubbles [U [2] containing a new phase. Unless the rate of bubble formation outpaces the expansion of the inflating 
r ■ background, the original phase is never completely consumed, and inflation becomes eternal. In this picture, many 
different phases can be seeded from an eternally inflating "parent" vacuum, leading to a patchwork-universe with 
diverse spacetime-dependent physical properties. For a review of eternal inflation, see e.g. Ref. |3j. 

Surprisingly, there may be directly observable consequences of living in an eternally-inflating multiverse. In eternal 
inflation, our observable universe resides inside one member of an ensemble of bubbles. These bubbles necessarily 
] collide, perturbing the homogeneity and isotropy of our own bubble interior, and providing a direct observational test 
^ of eternal inflation [3]. Assessing the likelihood of seeing bubble collisions, determining their effects on cosmology, and 
• • predicting their observational signatures has been the subject of a large body of work [5lI22j. This body of literature 
> has established the plausibility of models where collisions are likely, compatible with our observed cosmology, and 
leave observable signatures. However, an existence proof does not yet exist. For a review of much of this work, see 

Ref. Eg. 

Ci The primordial inhomogeneities left by bubble collisions are imprinted in the cosmic microwave background (CMB) 

radiation. The signature of an individual collision is an azimuthally symmetric modulation of the CMB temperature in 
a localized region. In general, the theory predicts a set of observable collisions, with the properties of each drawn from 
an in-principle calculable probability distribution. A search for the signatures of bubble collisions in the Wilkinson 
Microwave Anisotropy Probe ''2^ (WMAP) 7-year CMB data [25] was performed in Refs. [26l[27], which yielded an 
upper bound on the average number of observable collision signatures on the CMB sky. 

A key element in the analysis of Refs. [17] is the fact that each bubble collision is expected to leave a fairly 
generic set of signatures in the CMB jTH [151 HSl ES] , described by a phenomenological template with only a few 
free parameters. However, there is currently no way to connect the parameters in this phenomenological model 
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to the parameters in the fundamental theory. In order to go from the Lagrangian to the signatures of bubble 
collisions, one must determine the properties of the colliding bubbles, the immediate outcome of the collision, and the 
cosmological evolution of the perturbed bubble interior. Because of the non-linear nature of the field equations, a full 
and unambiguous determination of the cosmological signatures of bubble collisions requires numerical relativity. See 
e.g. Ref. 28 for a modern review of this field. 

In this paper, we present the results of fully-relativistic simulations of the collision between Coleman-de Luccia 
(CDL) bubbles containing a realistic inflationary cosmology. Although any given bubble undergoes many collisions, 
the future domains of influence of each of the collision events typically do not overlap until relatively late times inside 
the bubbles. It is therefore a good approximation to treat each collision as occurring in isolation, and we restrict 
ourselves to the study of collisions between two bubbles. 

The spacetime resulting from the collision between two CDL vacuum bubbles has S0(2,l) symmetry [TjiBJ, allowing 
us to treat the problem of bubble collisions in full generality using a I+ID simulation. This symmetry is expected to 
be mildly broken by fluctuations on the walls of the colliding bubbles [2UH?2j . However, expanding bubbles become 
more and more spherical (the perturbations remain much smaller than the overall size of the bubbles), and so our 
truncation to 1+1 dimensions should capture the most important dynamics. Einstein's equations in 1+lD are fully 
constrained by the scalar field configuration, allowing for a straightforward and efficient numerical implementation. 

The goal of the present paper is to study the outcome of bubble collisions and identify the properties that models 
must have in order to be consistent with our observed cosmology. This has been addressed to a certain extent in 
previous work using the Israel Junction Condition formalism J, 8j iTUl 1111 1131 122] , as well as in numerical simulations 
that neglected gravitational effects [TH [T71 [15] . In addition, Ref. [33] studied collisions in full numerical relativity for 
the collision of vacuum bubbles (see also Ref. [34] for a relativistic treatment of single bubbles) . Our fully relativistic 
treatment allows us to go beyond these previous studies and completely determine the phenomenology of bubble 
collisions given a particular model of the scalar field theory driving eternal infiation. In addition, we determine 
how well the approximations used in previous work hold up when compared against the fully relativistic solutions. 
In a forthcoming paper, we will extend this numerical framework to extract the observational signatures of bubble 
collisions. 

The remainder of the paper is organized as follows. We begin in Sec.|IT|by reviewing conclusions obtained from the 
Israel Junction Condition formalism and summarizing the set of questions we attempt to address with our simulations. 
In Sec. |III[ we perform the truncation to the 1+lD system of equations to be solved numerically, and discuss the 
generation of initial data. We then briefiy outline the numerical implementation and tests of the code in Sec. |V| and 
present the results of our simulations in Sec. |VI| More detailed treatments of a few technical aspects of the setup are 
presented in a set of appendices. We work in units where Mpi = 1 unless otherwise noted; Newton's constant is given 
by Gn = Mpi^. 

II. GENERAL PROPERTIES OF BUBBLE COLLISIONS 

In this section, we outline the general setup for bubble collisions using the Israel Junction Condition formalism [35j . 
This also serves to lay out a set of expectations for the numerics and illustrate the approximations made in previous 
work. After this introduction to bubble collisions, we summarize the main questions which are addressed by our 
numerical studies. 



A. The Israel Junction Condition formalism 

The generic bubble collision spacetime is shown in Fig. [Tj We work in the "Collision frame," the frame in which 
the two colliding bubbles nucleate simultaneously in some timeslicing | iTT| . It is always possible to perform a Lorentz 
boost to go to the collision frame in a spacetime with two bubbles. As shown in the figure, the two bubbles form 
(Regions II and IV in the figure) from the false vacuum (Region I in the figure) , expand, and then collide. The bubble 
interiors are perturbed in the future of the collision (Region III in the figure) ; the boundary of the affected region can 
carry energy density, and the interior of this region might be comprised of multiple subregions. We refer to the bubble 
on the left, which we assume to contain an inflationary cosmology consistent with our observable universe, as the 
"Observation bubble." The "Collision bubble," which we always place on the right, can have more general properties, 
and serves to perturb the Observation bubble interior. 

In general, there are a number of possibilities for the region to the future of the collision (Region III), and the 
interfaces separating this region from each bubble interior: 

• Hill = Hiy: Region III belongs to the Observation bubble; the c interface is a shell of radiation required to 
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FIG. 1: The generic bubble collision spacetime in the thin- wall approximation is split into four regions. Region I is the false 
vacuum, Region II is the interior of the Collision bubble, Region III is the future domain of influence of the collision (which 
might be composed of further sub-regions), and Region IV is the interior of the Observation bubble. In the Israel junction 
condition formalism, the Observation and Collision bubbles are separated from the false vacuum by infinitesimally thin walls 
(a and d) characterized only by their tensions. In addition, Region III is enclosed (depending on the context) by either domain 
walls or shells of scalar radiation (b and c), each characterized by their tension and energy density. 



conserve energy and momentum; the b interface is a domain wall separating the interior of the Observation and 
Collision bubbles if they are different vacua, or a shell of radiation if they are identical vacua. 

• Hill — Hii. Region III belongs to the Collision bubble; c is a domain wall if the Observation and Collision 
bubbles contain different vacua (it is a shell of radiation if they are identical); 6 is a shell of radiation. 

• Hill 7^ Hii, Hi\: Region III is in a distinct vacuum, different from that of the Observation or Collision bubbles; 
h and c are domain walls separating Region III from the Collision and Observation bubble interiors. This is 
known as a classical transition [T71 [TS] . 

The collision spacetime in the Israel junction condition formalism is constructed by matching vacuum solutions of 
Einstein's equations across interfaces characterized by some tension Ui (where i — a,b,c,d denotes the walls separating 
the various regions in Fig. [T]). Each region corresponds to a different vacuum of the underlying potential, and the 
tensions are determined by the potential barriers between vacua. As discussed in the introduction, the collision 
spacetime has S0(2, 1) (hyperbolic) symmetry. Using a hyperbolic version of Birkhoff's theorem [51 [TOl EJ [T3], we 
can find the most general vacuum solutions to Einstein's equations possessing this symmetry: 

ds^ = -A^{z)-^dz^ + A^{z)dx'^ + z'^dHl , (1) 

where = dx^-l-sinh^ xc?</'^ is the metric on a unit hyperboloid, and the coordinates generally range from < z < cx), 
— oo < X < oo, < X < oo, and < < 27r. The metric functions A^{z) are 

A,{z) = l-'^-^+H^^z\ (2) 

where we have included the subscript 7 = 1, II, III, IV to label the different regions in the collision spacetime. There 
are two constants in each metric function, the mass parameter (akin to the Schwarzschild mass parameter) and 
the Hubble parameter H^ (related to the cosmological constant by A = 3H^). When — 0, the metric describes 
a hyperbolic foliation of de Sitter space (hyperbolic de Sitter, which we denote by HdS). In our conventions, H^ > 
(< 0) for positive (negative) vacuum energy. In the following, we only consider the case where H^ > 0. 

Returning to Fig. [T] we can identify this sketch as a slice through the full bubble collision spacetime in the x — z 
plane defined by x = and 4> = const. All the plots shown in later sections are along a similar slice through the 
collision spacetime. Integrating Einstein's equations across each interface and requiring continuity of the metric yields 
an equation of motion for each interface. The continuity of the metric across each interface requires z to be continuous. 
However, the x coordinate is not continuous, and we therefore have to take into account different coordinates, xl and 
Xr, on either side of each interface. Defining ki = 47rcri, the equations of motion in terms of the proper time measured 
by observers along the interface are given by: 
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where and An are the metric functions on the left (L) and right (R) side of a particular interface in Fig. [Tjand 

^« - 2^, ■ 

The sign of xl,r determines the direction in which the wall is moving from the perspective of observers on either side. 
From Eq. [2] for the metric functions, there are two possible behaviors for the wall at large z: normal and repulsive. 
Normal walls require — > k^, in which case i has the same sign on both sides of the wall, and the wall 
accelerates towards the side with the larger Hubble parameter. Repulsive walls occur for — < k^. In this 
case, i takes opposite signs on either side of the interface, implying that it is accelerating away from both sides. 
This behavior arises due to the gravitationally repulsive nature of domain walls [361 137| . The simplest application of 
these equations of motion is to calculate the trajectory of the wall (of the normal type) of a single bubble from the 
perspective of the false vacuum. This is given by 

z = Hp^ [(1 - HIrI) sec^ (Hp [x + Ax])] , (6) 

where Hp is the Hubble parameter in the false vacuum from which the bubble is formed, and Rq is the critical radius 
of the bubble in the thin wall approximation. More generally, Eq. [3] determines the motion of a post-collision domain 
wall, should one be produced by the collision. 

The kinematics of a collision are fully determined by the location Zc of the collision. This is equivalent to specifying 
the distance between the bubble centers in the Collision frame. The surface tensions of the interfaces and vacuum 
energies are determined by the underlying potential, but the mass parameters in Eq. [l] for each region remain un- 
determined. Imposing the condition that there are no conical singularities at Zc (equivalent to imposing energy and 
momentum conservation [38^) yields one equation relating these mass parameters to the vacuum energies and tensions. 
In the limit where the incoming and outgoing interfaces are highly relativistic, this relation takes a particularly simple 
form |S] , which with four distinct post-collision regions as shown in Fig. [l] is given by 

Ai{z^)Aui{zc) = Au{zc)Aiv{z^). (7) 

For vacuum bubbles, the mass parameters in regions I, II, and IV are zero. Assuming that region HI is not subdivided, 
there is only one mass parameter, which is determined by Eq. [7] 

z'c Hi, - gFv + H^i (1 + Hl,zl) Hi (1 + gfVzg) 

= y i^Hfzi • 

Although the mass parameter increases without bound as the collision occurs at increasing Zc, the collision between 
two bubbles can never produce a black hole [7| [HI HO] ; this would be indicated if there was a solution with ^ni = 
which, from Eq. [Tj is not possible. 

The thin- wall solutions provide an excellent way of characterizing the possible outcomes of the collision between two 
vacuum bubbles. However, these solutions do not take into account field dynamics, and therefore are not sufficient to 
compute the effects of a bubble collision on the cosmological evolution inside a bubble. Before tackling this question 
using numerical methods, let us first summarize some of the specific questions we address. 



B. Questions for the relativistic simulations 



The set of questions we are primarily interested in answering using our numerics are as follows: 

• How valid is the thin-wall approximation, and does the Israel junction condition formalism miss any dynamics 
that are important for determining the outcome of bubble collisions? How important are the internal dynamics 
of the walls themselves? 

• In the thin-wall approximation, the interior of the colliding bubbles is treated as having constant energy density. 
Since one would ultimately like to embed a cosmology into the bubble interiors (to make contact with our own 
observable universe), this is clearly a key limitation. Previous work used the thin- wall solutions as a background 
in which to evolve a scalar field fH' but to date, there has been no full solution of cosmological evolution inside 
colliding bubbles. Taking into account full General Relativity, how does the collision affect cosmology? 
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• Non-gravitational simulations predict that the outcome of a bubble collision might depend sensitively on 
the underlying potential. In particular, Ref. pjj found that models with inflaton potentials of the "small-field" 
type (those driving inflation near a maximum or inflection point in the potential) were far more sensitive to the 
effects of bubble collisions than models with inflaton potentials of the "large-field" type (such as inflation in a 
quadratic potential). If true, this would have important implications for the connection between the possible 
presence of bubble collisions and other cosmological observables, in particular the tensor-to-scalar ratio [39j . 
Can inflation of the small-field type occur to the future of a bubble collision? What properties must a potential 
have in order to produce collisions that are consistent with our observed cosmology, yet produce observable 
signatures? 

• Depending on the structure of the potential underlying eternal inflation, a collision between two bubbles can 
push the field into the basin of attraction of a new vacuum, causing a "classical transition" to a new phase [IT] . 
In the absence of gravity, classical transitions can only give rise to lasting regions with a lower energy density 
than the colliding bubble interiors. However, in the presence of gravity, the thin-wall analysis of Ref. [52] showed 
that it is kinematically possible to create lasting regions of higher energy density in the collision between two 
vacuum bubbles. Such transitions can even prevent a portion of the interior of a bubble undergoing a big crunch 
from encountering a singularity. Do these result hold true once the detailed field dynamics are included? 



III. EQUATIONS OF MOTION 



We now consider non-vacuum solutions of Einstein's equations, and derive the equations of motion used in our 
implementation within the Cauchy approach (see e.g. Ref. [28j). Adopting coordinates = {z,x^) — {z,x,x,<P), the 
most general line element with hyperbolic symmetry can be written as: 



{-a^ + a^j3'^)dz^ + 2a'^(3dzdx + a^dx"^ -f z^lF' [dx^ + sinh^ 



(9) 



with {a^,6^} metric variables describing the intrinsic curvature 7ij = diag[o^, sinh^ x] ol (spacelike) hyper- 

surfaces defined by z = const and {a,/3} coordinate conditions. It is convenient to introduce the extrinsic curvature 
—2aKij = {dt — to express the Einstein equations within the ADM formalism [JD]. The equation governing 

the behavior of the scalar field is 
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We reduce this equation to first order form by introducing the variables 



n = - {ip-Pv') 

a 



(10) 

(11) 
(12) 



The general equations of motion for hyperbolic symmetry and arbitrary {a,/3} are derived in Appendix [A| For our 
current purposes we exploit the coordinate freedom to set 



1, /3 = 0, 



(13) 



motivated by the vacuum solutions to Einstein's equations with metric Eq. [T] With this choice the Einstein equations 
(of motion and constraints) result in: 



a 
b 



z 



IGtt 



2a2 
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1 fa' 
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8ttV 
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The remaining freedom can be tied to the choice of the lapse function. Again, a particularly simple option is to define 

^ (20) 



which, when substituted into the momentum constraint equation (Eq. 19 1, gives 



aM2j 



(21) 



where we have restored factors of Mp\. 



Next, we can solve for K'^ x in the Hamiltonian constraint (Eq. 18), 



(22) 



and replace it in the time-evolution equation for a to obtain (restoring factors of Mpi) 



2z 



2a2 



V{ip) 



(23) 



Now, an evolution equation for a can be defined by inserting K^^^ and into the time-evolution equation for K^^. 
Restoring factors of Mp\, this is given by 

If ^. Sn 



2z 



a~ \ 1 + I 



27rza^ 



In addition, we have the equations of motion for $, 



$ = a 



(24) 



(25) 



and n. 



n 



-n 



a a a\ a , 

Tr]<^+ -aad^V . 

a a-^ a 



(26) 



For the simulation, it is convenient to work with a set of units where time, distance, and energy are measured in 
terms of the false- vacuum Hubble constant Hp , defined by 



„2 _ StT 



We now define the following dimensionless variables: 



(27) 



ip ~ $ . n ~ V 

z = Hfz, x^Hpx, — , $ ^ , n= , V ^ 

Mpi Mp\Hp Mp\Hf M^Mi 



(28) 



The evolution equations become: 



da a 
dz 2z 



1 - 1 + Svrz^ 



n2 + $^ 

2a2 



(29) 



da a r 
25 



1 



{l+STlZ^V{(p)^ C? 





^T^ZC? 

+ 2- 













(30) 



dn 

dz 



z 



1 da a da \ ^ a d$ dy 

a dz a? dz J a dz dip ' 



(31) 
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dz ( a ~ ) ^ ^ 



In the following, we suppress the tildes for notational clarity. 



The system of equations (29 31) fully determine the future evolution of the field and metric variables given consistent 
initial data (e.g. satisfying the constraint equation Eq. [21] ). The simplest test of the validity of these equations is to 
reproduce the vacuum solutions (Eq. [I]) in the limit where all derivatives of the field are zero and V = const. We 
demonstrate this in Appendix [BJ It is also illustrative to consider the evolution of fluctuations in the field in the false 
vacuum background. This is outlined in Appendix [C] where we show that the characteristic speed of fluctuations is 
given by: 

c = ±(l + i7|z^)"\ (33) 
which ranges from unity at small z to zero as z becomes large compared to Hp. 



IV. DETERMINING THE INITIAL DATA 



Our strategy for defining the initial data for bubble collisions proceeds in two steps. First, given a potential for the 
scalar field, we construct the CDL instanton for each possible transition from the false vacuum. A slice through the 
CDL instanton fully determines the initial data for a single-bubble spacetime. When the bubbles are small compared 
with the horizon size defined by the false vacuum, and have tensions small compared with the energy scale of the 
false vacuum, they introduce small perturbations to the false vacuum background. In this limit, we can construct the 
two-bubble spacetimes by linear superposition. Next, we need to express the field and metric configuration defined by 
the CDL instanton in terms of the variables used in the simulation. To do so, we adopt an initial hypersurface at z = 0, 
where the a;-coordinate in the hyperbolic slicing can be identified with a proper distance. It is then straightforward 
to map the field configuration from the CDL instanton onto the z = hypersurface, and perturbatively determine 
the field and metric functions at z = dz. This provides consistent initial data to evolve the solution forward from 
z = dz: we conflrm the validity of our initial data by evaluating the Hamiltonian and momentum constraints and 
checking for the convergence of the full solutions. In what follows, we briefly describe the procedure adopted and 
defer to Appendix [P] for further details. 



A. Solving for the CDL instanton 

The CDL instanton [TJ [5] is a compact solution to the Euclidean equations of motion with Euclidean metric: 

ds^ = dt^ + p{tfdnl . (34) 

The metric function p{t) has two zeros at t = and t — tmaxi and evolves according to the Euclidean Friedmann 
equation 

As described in Appendix [P] (and further in Refs. [ITl I42j). the importance of gravity in the computation of the 
instanton is characterized by the parameter 



SttM^ 



where M is a scale parameterizing the width of the potential barrier. In the numerical solutions that follow, a typical 
value is e = 0.01. With the metric Eq. [34j the field equation is: 

d^ idp^dip ^ dy_ 

df^ p dt dt dip' ^ ^ 

The CDL instanton is a non-singular solution to the equations of motion that interpolates between the basins of 
attraction of the true and false vacua. Non-singular solutions satisfy the boundary conditions: 

(^(t = 0) ~ </JT, 5(t = 0)=0, p(i = 0)=0, ^{t^t^^^)^^F, $(t = tmax)=0, p(i = t^,,) = 0. (38) 

dt dt 
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FIG. 2: The numerically generated instantons for a potential of the form Eg. |D6| with a = 0.1 for e = 0.0075. The dimensionless 
variable x is related to the field Lp, r to the metric function p, and s to f as defined in Eq. |D2[ The prediction for the critical 
radius (based on Eq. D20 in the thin-wall approximation) is shown as the blue dashed line. 



A review of CDL instantons and their construction can be found in e.g. Ref. [43] . 

In general, it is not possible to solve analytically for the CDL instanton (although see Ref. [H] for some special 
cases where an analytic solution can be found). We therefore adopt a numerical scheme to solve the double-boundary 
condition problem defined by Eq. |38[ We describe our numerical methods in Appendix [Dj An example output is 
shown in Fig. [2j 



B. Obtaining initial data in the hyperbolic slicing 

We are now in a position to determine the initial data for the field and metric functions in the hyperbolic slicing used 
in the simulation. We first discuss how initial data for a single bubble is constructed, and then add a colliding bubble 
to the picture. The CDL instanton can be analytically continued to obtain the post-tunnelling field configuration 
^p{t). The metric is given by: 

ds^ = dt^ + p{tf + cosh^ mdVl?^ . (39) 

In the absence of a bubble, the slice at 5* = corresponds to the throat of de Sitter. When a bubble is present, for 
small e, this slice is a small perturbation to the throat of de Sitter space (see Ref. [23] for more detailed discussion of 
the coordinates relevant to the CDL bubble spacetime). Here, t measures a proper distance. In the hyperbolic slicing, 
z = can also be identified with a portion of the perturbed throat of de Sitter. In the limit where z — 0, the metric 
is given by 

ds^ = -dz^ + dx^ + z^dHl . (40) 

Therefore, since both x and t measure a proper distance on the throat of de Sitter, we can simply substitute x for 
t in the instanton solution. This is the field configuration on the z = slice. Further, since the curvature of the 
transverse components of the slices dominates the energy density at z = 0, the metric functions do not concern us 
here. Examining the equations of motion, there are potentially various singularities at z = 0. We must therefore find 
the initial conditions after one time step dz. To do so, we first consider the necessary asymptotics for the various 
functions. The field ip is expanded as 

(f ^ ifio{x) + ip2{x)z^ , (41) 
where the first order term is excluded since we require dip/dz = for a non-singular solution. We also have 

a = 1 + a2{x)z'^ , a = 1 — a2{x)z'^ . (42) 
Using these asymptotic relations, and working to lowest order in z, we have: 



(43) 
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Now, substituting into the equations of motion (Eqs. 29-32), we can solve for the unknowns (p2{x), 02 (a^), c(2ix), 
in terms of ipo^x). From the constraint equation for a, 



da2 2 a dipo 2 
ax ax 



while from the evolution equation for a, 



2zq;2 



a2 - ^ttV + 2tt 



dx 



3 3 \ dx 



This matches with the expectation for a constant potential V in the small-z limit: 

4:TTV 



"HdS = 1 - 



-z' + .. 



From the evolution equation for a, we have 



2a2Z 



AttV - q;2 + 27r 



dm 
dx 



We also have the field equations. From the evolution equation for H, we find: 



2^2 



dx"^ 



Atp2 - dV, 



dx"^ 



(44) 



(45) 



(46) 



(47) 



(48) 



The evolution equation for $ yields an identity. Note that substituting the definition of ip2 into the constraint equation 
for a (Eq. 21 ), we find consistency. Substituting the expressions for ip2{x), a2{x), a2{x), and 1^30(2;) into Eqs. 42 and 43 



we set initial data at z = dz <^ 1. 

The above procedure is perfectly rigorous for initial data involving a single bubble. In order to discuss configurations 
with two bubbles, we must make a few approximations. Assuming that the bubbles represent minor perturbations to 
the false vacuum background, we determine the initial data for two different bubbles independently, add the two field 
configurations, and calculate the quantities in Eqs . [43| and [42| as above. To the extent that each of the bubbles evolves 
before the collision just as it would in isolation, this is a good approximation. In practice, this amounts to choosing 
potentials where one instanton endpoint is very close to the false vacuum, the initial radii of the bubbles are small 
compared to Hp^, and ensuring that the bubbles are sufficiently far apart on the z = hypersurface for the field to 
be very close to the false vacuum between the bubbles. For bubbles satisfying the thin-wall approximation described 
in Appendix [D] these conditions are almost always satisfied. Our choice of e ^ 1 ensures that we are well within this 
regime. 



C. Potentials 



For our simulations, we must supply a set of scalar potentials. We define a class of piecewise potentials constructed 
in a manner similar to the procedure outlined in Ref . |12j . There are 4 possible segments for each potential, as shown 
in Fig. [sj These segments correspond to two potential barriers (labeled Ti and T2 in the figure) and two slopes leading 
away from each potential barrier towards a set of minima (labeled Ci and C2 in the figure). 

By convention, T2 and C2 correspond to the portions of the potential relevant for the Observation bubble. For 
studies of the collision between two vacuum bubbles, we take Ci and C2 simply to be minima of the potential. For 
bubbles with a cosmology, the C2 portion of the potential is chosen so that an undisturbed Observation bubble 
contains an epoch of inflation which lasts at least 60 e-folds. We consider two broad classes of inflationary potentials 
C2: those for which inflation occurs over a field range larger than or comparable to Mpi (large-field inflation) and 
those for which inflation occurs over a fleld range smaller than Mpi (small-field infiation). 

The functional form for the various segments is given by 



V 



Vci ii (f < fTi 

Vti if y^Ti < <(5 < 

Vt2 if < (/3 < ipT2 

Vci if (PT2 < V 



(49) 



10 



V 




FIG. 3: The generic form of the potentials used in our simulations consists of four components. Tl and T2 describe the 
potential barriers between the three vacua of the potential. CI and C2 describe regions of the potential leading to the vacua 
and can drive a period of inflation. 



where, when there is a cosmology inside the Collision bubble, we have 



In all cases, we have 



Vr.^-'^i^- + - ^.f + + M^C, , (51) 



Vr.^-^i^~ ^.f - - ^.f + + M^C, . (52) 



For large-field inflation, we assume a form of the potential given by: 



2 



2 



while for small-field inflation, we use: 



Vc2 = Vt2 - - ^T2f + ^ (<y3 - ipT2)^ ■ (54) 

Here, the parameters we must specify are ai, a2, M, Vqi, ipmi, ^m2, Vjii ^j2, Ci (and for the small-field model, 
V and A). This corresponds to specifying the false vacuum energy density, barrier widths, and the length of the 
infiationary regions of the potential on either side. All other parameters are fixed by requiring continuity of the 
potential and its derivatives. We have simulated a wide variety of potentials, but consider only a representative 
sample here. The parameters defining each of these examples are recorded in Appendix [E] The potentials giving rise 
to vacuum bubbles are shown in Fig. |4j the potentials giving rise to bubbles with an interior cosmology are shown in 
Fig. [5] These particular choices are motivated in the following sections. 



V. NUMERICAL IMPLEMENTATION 



A. Testing for convergence 



In this section, we discuss the details of our numerical solution of the equations of motion (Eqs. 29 32). Since the 
bubbles we consider are small perturbations to the false vacuum, we can use the behavior of null geodesies in HdS 
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FIG. 4: Potentials giving rise to vacuum bubbles. For the potentials on the top row, the false vacuum is located at (p = 0, and 
there are two possible CDL solutions, connecting the false vacuum to positive or negative ip. In both cases, the lowest energy 
vacuum is the one at positive p. The potentials on the bottom row have a false vacuum at negative </?. There is only one CDL 
solution, connecting the false vacuum to the vacuum at ip — 0. The lowest energy vacuum for the potential V3 is at positive p 
while the lowest energy vacuum for the potential V4 is at = 0. 



to determine the spatial size of the simulation. This is discussed in Appendix [Fj where we show that a box of size 
X € [— 7r/2, Xc + 7r/2], with Xc denoting the center of the CoUision bubble, contains the two colliding bubbles in their 
entirety. For larger bubble separations, we also consider a box extending from the origin of the Observation bubble 
at a; = to the origin of the Collision bubble at Xc- 

To integrate the equations of motion forward in time, we use the method of lines on a fixed uniform grid with a 
4th order Runge-Kutta scheme (see, e.g. Ref. [IS]). To discretize the equations in space, we use Finite Difference 
approximations for the spatial derivatives (satisfying summation by parts) which are 4th order accurate at all inte- 
rior points, with the exception of two points at the boundaries where the convergence is 2nd order. The resulting 
implementation, as discussed in Refs. |451 146] . ensures a consistent treatment of the problem and guarantees stability 
in the linearized case. The obtained solutions should converge to 3rd order accuracy in the discretization length. We 
confirm this through several tests, which we now describe. 

We choose the potential V2 shown in Fig. |4j and set up initial conditions consisting of two different bubbles with 
the Collision bubble located at HpXc — 0.75. A contour plot of the field (p is shown in Fig. [6j The details of the 
solution are not immediately relevant, but this is a fairly representative example of a collision spacetime. We monitor 
the solution up to z = 4.3, which includes the collision and the post-collision domain wall. We evolve this scenario on 
a uniform grid with number of points Np = 2^ x 10'^ (p = 0..5) and compare the obtained solutions. In particular 
we monitor the (^2 norm of the) difference of solutions with successive values of p (i.e. for a given field, we define 
Cp-p+i = \\ipp — (pp+i\\2)- These are shown in Fig. [t] (left panel) highlighting how, as resolution is increased, the 
difference decreases as expected for a convergent solution. We can further test the convergence rate q as estimated by 

29 ^ Wfp " y^P+ilb ^ ep-p+i ^^^^ 
ll'/'p+i - 'PP+2II2 ep+i-p+2 

The results are shown in Fig. [7] for the field ip, indicating that the solution converges; further, as the resolution is 
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FIG. 5: Potentials giving rise to bubbles with an interior cosmology. The potentials on the top row are of the large-field type, 
both have the false vacuum at = 0, and both allow for two possible CDL solutions connecting the false vacuum to regions of 
positive or negative ip. The potential LI has cosmological evolution inside both bubbles, with a much steeper slope at negative 
ip than positive ip. We show both the full potential, and a plot zoomed in around the position of the false vacuum. The potential 
L2 omits the cosmological evolution at negative ip, but has the same potential at positive ip as LI. We show only the portion 
of the potential L2 in the vicinity of the false vacuum. The potentials on the bottom row are of the small-field type, and again, 
have the false vacuum at <^ = 0, with two possible CDL solutions. The potential SI is defined using Eq. [54]while S2 is defined 
using Eq. |53| inflation occurs in the vicinity of an inflection point for SI while it occurs in the vicinity of a hilltop for S2. 

increased, the overall convergence rate tends to the expected value of 3 [46 . In addition, we have confirmed that the 
constraint Eq. [2T]is wcW preserved throughout the simulations. The resolutions yjc use here are representative of what 
is computationally feasible with our current implementation (the run-time for the highest resolution shown here was 
approximately 48 hours on a single core). 

Unless otherwise noted, all of the simulations in the following sections are run with a total of = 32 x 10'^ spatial 
grid points. Having demonstrated a convergent code, we now perform a number of cross-checks of the expectation 
from the thin-wall single bubble-solutions outlined in Sec. |II A[ 
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FIG. 6: A contour plot of the field tp[x, z) for the collision of the two different types of bubbles allowed by the potential V2 
shown in Fig. [4] The color bar indicates the value of the field at each point in spacetime. In this case, the post-collision domain 
wall is repulsive, and accelerates away from the interior of both bubbles. 
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FIG. 7: Demonstrations of the convergence of our code for the simulation shown in Fig. [6] On the left, we show the absolute 
value of the difference between solutions with discretizations obtained for p and p+ 1, illustrating how the difference decreases 
with increasing resolution as expected for a convergent implementation. On the right, we plot the estimated convergence rate 
as calculated from Eq. [55] as the resolution increases, the value obtained approaches 3, which is consistent with the expected 
value for our implementation. 



B. Simulation of a single bubble 

In this section, we study the evolution of single bubbles. Our purpose is threefold: we check the numerical solutions 
against the predictions of the Israel junction condition formalism outlined in Sec. |II A[ discuss the general features of 
the coordinates used in our simulation; and illustrate the limitations of our numerical implementation. 

The expectation from the Israel junction condition formalism is that the walls of vacuum bubbles follow the 
trajectory given by Eq. [6j The left panel of Fig. [8] shows a contour plot of the field for a single Observation bubble 
generated from the potential VI. The field is in the vacuum at </? = everywhere outside the bubble (dark, purple 
area), and in the vacuum at positive everywhere inside the bubble (light, green area). Even though the bubble wall 
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FIG. 8: The numerical simulation of a single Observation bubble formed from the potential VI. On the left, we show a contour 
plot of ip{x,z). As this is a vacuum bubble, the field is constant on either side of the bubble wall. The yellow dashed line is 
the wall trajectory predicted by the thin-wall approximation, illustrating the excellent agreement between the numerical and 
junction condition solutions. In the middle, we show a contour plot of the numerical solution for a{x, z) (black solid lines) as 
well as the junction condition solution (red-dashed lines). The agreement is excellent, except in the vicinity of the bubble walls 
at large z. On the right, we show a contour plot of the numerical solution for a{x, z) as well as the junction condition solution 
(red-dashed lines). The disagreement between these two solutions is due to the fact that in the junction condition formalism, 
the x-coordinate is discontinuous across the wall, while in the simulation it is treated as a continuous coordinate. The solutions 
agree in the bubble center and the bubble exterior, which is expected. In the numerical solution, a must continuously interpolate 
between the solution in the bubble center and the bubble exterior. 



is expanding at an increasing velocity as time goes on, it reaches a constant comoving size. Since x is a comoving 
coordinate, the bubble walls asymptote to a constant position as z — >■ cxd (as can be seen from Eq. p|. This can also 



be understood from the fact that the characteristic speed of fluctuations given in Eq. 33 goes to zero as z — >■ oo. 
In Fig. [8] the thin-wall solution is over-plotted as the dashed yellow line, showing that the agreement between the 
numerical and thin- wall solutions at this resolution is excellent. 

We can also determine how the metric functions a and a compare with the expectation from the junction condition 
solutions. Identifying a = A-^ from Eq. [ij the thin-wall solution has a = (1 + Hfz'^)^'^ outside the bubble and 
a = {1 + Hf^z^y^^ inside it (using the conventions of Fig. [T]). In the middle panel of Fig.|8]we show a contour plot of 
a (solid black lines) along with the junction condition solution (red-dashed lines). The agreement is good, except at 
large z in the vicinity of the wall. We return to this below. In the case of the metric function a, the comparison is not 
as straightforward. In the junction condition solution, the x coordinate is discontinuous across the bubble wall, while 
for the numerics we have (appropriately!) chosen a coordinate system in which x is continuous. A contour plot of a 
is shown in the right panel of Fig. m (solid black lines), and compared to the expectation from the junction condition 

— 1/2 

solution (red-dashed lines), where a = Aj . Outside the bubble, and at its center, the thin- wall and numerical 
solutions match, as they should. Between the bubble center and the wall, the behavior of a differs from the junction 
condition solution in order to enforce the continuity of x across the wall. 

We now turn to bubbles containing a cosmology. As a representative example, we simulate a single Observation 
bubble in the potential LI. A contour plot of the field is shown in the left panel of Fig.[9j The metric functions a and 
a are qualitatively similar to those shown in Fig. [Sj since in this example the field is undergoing slow-roll inflation 
inside the bubble, and the solution (over the timescale of the simulation) is therefore quite close to pure HdS. The 
theoretical expectation for the evolution of the field inside the bubble is found by integrating the field and Fricdmann 
equations for a scalar field in an open FRW universe, with an initial condition set by the CDL instanton. In our 
numerical solution, the field at the center of the bubble as a function of the elapsed proper time t = J dz a{z, a; = 0) 
should match this theoretical expectation for the time evolution of the field. This is indeed the case, as shown in the 
right panel of Fig. [9j Moving away from the bubble center, we should find that the surfaces of constant field extend 
all the way up to z ^ cx). This would be true if the bubble interior truly contained an open FRW universe. The 
numerical solution shown in the left panel of Fig. |9] does not have this behavior at large z, although it converges to 
it with increasing resolution. 

In our simulations of both vacuum bubbles and bubbles containing a cosmology, the numerical solutions are worst 
near the bubble walls. Physically, this is to be expected. As the bubble expands, the wall becomes thinner and 
thinner due to Lorentz contraction (the velocity of the wall is continually increasing) and the fact that a; is a comoving 
coordinate. In Appendix |Fj we show that if the metric is approximately HdS, the width of the wall in terms of the 
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FIG. 9: Numerical simulation of a single Observation bubble formed from the potential LI. On the left, we show a contour 
plot of ip{x,z). The trajectory for the wall in the junction-condition solution (yellow-dashed line) can be found under the 
assumption that the bubble interior has a constant vacuum energy; the agreement with the numerical solution is excellent. 
The field is in the false vacuum outside the bubble, and evolves on the bubble interior down the potential towards increasingly 
positive values of ip. The numerical solution has the correct qualitative behavior except in the vicinity of the bubble walls. 
As the resolution is increased, we expect that the contours in the bubble interior to go up to z — >■ oo. In the right panel, we 
plot the field in the center of the simulated bubble ip{x = 0, z) as a function of the elapsed proper time (black-solid line). The 
slow-roll solution obtained by integrating the field and Friedmann equations is overplotted (dashed-yellow line); the agreement 
is excellent. 



simulation coordinate x is inversely proportional to z. For our standard resolution of 32 x 10'^ spatial grid points, 
and for the potentials we study (which fixes the initial width of the bubble walls) , this imposes an upper limit on the 
run-time of roughly z 10 for our simulations. In fact, our ability to accurately track the solution in the vicinity 
of the walls breaks down somewhere in the range between z — 2 — 4 at the standard resolution, depending on the 
potential (see the contour plots in Figs. [8]and|9|. Thus, in order to accurately model bubble collisions, we must 
require that the collision occurs before our description of the walls breaks down, and we always impose Zc < 4. 

The numerical solutions we obtain still exhibit an overall convergent behavior, and so Richardson extrapolation |47j 
could be employed for better accuracy. Nevertheless, it is clear that for even later times the solution drifts away 
from the convergent regime, which limits the application of this approach. In addition, computation time restricts 
the number of grid points to perhaps two to four times our standard resolution, resulting in an increase by only a 
factor of two in run time. Thus, the main limitation of our current implementation is the restriction on run time, 
and the associated restriction on the kinematics of the collision. Because significantly longer run times are necessary 
to accurately extract the cosmological signatures of bubble collisions, we restrict ourselves in the current paper to 
exploring the immediate outcome of bubble collisions. To move beyond this, we are extending our implementation 
to exploit adaptive mesh refinement and alternative coordinate conditions (for lapse and shift). Our results on the 
cosmological signatures of bubble collisions will appear in a forthcoming publication. 

VI. RESULTS 

We now discuss the results of our simulations, beginning with the collision between vacuum bubbles (i.e. those not 
containing the inflationary plateaus CI and C2 in Fig. |3]), and then considering the outcome of collisions between 
bubbles that contain an interior cosmology. 

A. Colliding identical vacuum bubbles 

In this section, we discuss the collision of identical vacuum bubbles generated from the potentials shown in Fig. |4] 
Since the vacuum in each of the colliding bubbles is identical, the two bubble interiors must eventually merge. 
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FIG. 10: A contour plot of the field for a collision between identical bubbles formed in the potential VI. The solid lines and 
colorbar indicate the value of the field. In the immediate aftermath of the collision, a pocket of the false vacuum is formed. 
This pocket expands, re-collapses, and then re-forms. After a few cycles of this, a coherent, oscillating field configuration known 
as an oscillon is formed. 



However, to satisfy energy-momentum conservation, the energy of the collision must be dissipated in some way. 
The form of energy loss is an input to the junction condition formalism described in Sec. |II A[ and it is typically 
assumed [71 [HI [IHl [HI [H] that the energy of the collision is released in the form of scalar radiation. However, there 
are many other possibilities once the full dynamics of the field and metric are considered. 

In Fig. [TOj we show a contour plot of the field for the collision between two identical vacuum bubbles formed 
from the potential VI with an initial separation of HpXc = 1. Contour plots of the metric functions are shown in 



Fig. 11 In the immediate aftermath of the collision, a small pocket containing the false vacuum is formed. This 
region expands, re-collapses, and forms again when the walls re-collide ^ . In addition, some (relatively small) amount 
of scalar radiation is emitted at each collision. After a few cycles of expansion and contraction, there is no longer 
enough energy for the field to return back up to the false vacuum. The energy does not dissipate, but instead goes 
into the formation of a coherent, oscillating structure centered about the location of the collision. Such structures 
are known as oscillons j48H53j . and arise due to non- linear interaction terms in the scalar potential. In this case, 
if we extend the solution along the symmetry direction, the oscillating field configuration extends indefinitely in x 
(the coordinate labelling the spatial position along the hyperbolae of constant z). The formation of oscillons in the 
collision between bubbles was first discussed in Ref. (see also Ref. [53 [5B]). For all sets of kinematics that we have 
studied, an oscillon ultimately results. This is quite different from the standard assumption, based on the thin-wall 
solutions, that the colliding domain walls eventually dissipate into scalar radiation. 

Let us now analyze each stage of evolution in this collision. In the immediate future of the collision, the expectation 
is that the "free-passage" approximation is satisfied [17L I18j . In this approximation, the bubble walls do not initially 
interact (as long as they are sufficiently relativistic) , but instead superpose linearly. Immediately after the collision, the 
field in a small region is at a position ip ~ 2(j)j'2 (twice the width of the potential barrier) . This is high enough on the 
potential that the field overshoots the minimum, bounces off the rising potential at positive ip, and forms an expanding 
pocket of the false vacuum. Because this pocket has higher energy than the bubble interior, it eventually collapses. 
If there is enough energy left, a new pocket of false vacuum is formed after the collapse of the first. This continues 
until there is no longer enough energy to create a pocket of false vacuum, at which point an oscillon forms. In Fig. |12[ 
we show the oscillon profiles on constant-z slices through simulated collisions with three different initial separations 



This behavior was first observed in Ref. fT", and is referred to in the Uterature on thin-waU sohitions as an Oscillatory solution (see 
Ref. [22]). 
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FIG. 11: A contour plot of the metric functions a (left) and a (right) for the collision shown in Fig. |10[ 
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FIG. 12: The field on constant z ~ 3 slices for simulated collisions between identical bubbles formed in the potential VI for 
three initial separations HpXc = 0.25 (black-solid line), HpXc = 0.50 (red-dashed line), and HpXc = 1.00 (blue dot-dashed 
line). The solutions have been translated so that the collisions each occur at Hpx — 0. On subsequent time-slices, the profile 
retains its shape, but oscillates in amplitude. Such structures are known as oscillons. 



{HpXc = 0.50, 0.75, 1.00). The mam effect of increasing the initial separation is a larger-amplitude, smaller-width 
profile. In each case, the profile retains its shape as z increases, but oscillates with decreasing amplitude. 

We can understand some quantitative aspects of the oscillon solution analytically. Neglecting the gravitational 
backreaction of the field on the metric, Eq. |C1| gives the field equation for (p. We assume that the solution can be 
written as a stably oscillating packet plus a correction: 



(p{z, x) = fo{z)go{x) + Sip{z, x) 



(56) 



Substituting this ansatz into the equation of motion, Taylor-expanding the potential about the true vacuum minimum, 
and keeping only the lowest-order terms, we have: 



dzzfo - 

In the limit where z 3> H^^, 
the equation of motion is: 
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(57) 



we can neglect the contribution of the gradients in the field, and the general solution to 
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(58) 



Therefore, the expectation is that we have fluctuations of decreasing amplitude and increasing period as a function of 
z. The 5lp{z, x) term in our ansatz captures any additional effects, such as the dissipation of the oscillon into scalar 
radiation. 
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|12| we plot the field as a function of z. 
The qualitative agreement is quite good, although the analytic 
solution does not get the drift of the frequency of oscillation exactly correct. As the initial bubble separation is increased, the 
packet amplitude becomes larger. 



FIG. 13: Cutting through the center of the intermediate o scil lon profile shown in Fig. 
Overplotted (red-dashed lines) is the analytic solution Eq. 



Cutting through the oscillon profile at constant x in our simulations, we show in Fig.fTslhow the amplitude decreases 



with time for the intermediate case plotted in Fig. 12 We also plot the solution Eq. 58 (red dotted line) obtained in the 
limit where the gravitational back-reaction of the field configuration is negligible, matching the analytic and numerical 
solutions at Hpz = 2. First, note that the analytic solution approximates the z-dependence of the amplitude of the 
oscillations very well. This indicates that the dominant form of energy-loss for the oscillon is due to the background 
expansion, and not scalar radiation being shed from the configuration. The frequency of the oscillations is initially 
matched quite well by the analytic solutions, but there is a mismatch at larger values of z. This could be due to 
gravitational effects not captured by our analytic approximation, or non-linear corrections to the period of oscillations 
about the potential minimum. Nevertheless, the qualitative agreement between the analytic and numerical solutions 
is quite good. We also study the field as a function of proper distance on slices at z = 5 and z — 10, observing that the 
oscillon profile roughly doubles in width over this timescale. Such behavior is consistent with the profile undergoing 
expansion with the background, since the physical distance between two points of constant x has also doubled. The 
oscillons formed in these collisions are therefore not gravitationally bound. 

In the thin-wall approximation, the mass parameter in the metric function Eq. [2] characterizes the effect of the 
collision on the metric describing the bubble interiors. Referring back to Fig. [l] we can identify Region III with the 
false vacuum, use Eq. |8]to predict the mass parameter, and directly compare the metric function a between the thin- 
wall and numerical solutions. In the numerical solution, at the location of the collision, the metric function should 



jump, and then be dominated by the 1/z behavior dictated by the mass parameter. In Fig. 14 we show a at the 
location of a collision between two identical vacuum bubbles with an initial separation HpXc = 0.75 (left panel), and 
a at a location which intersects the outgoing shell of scalar radiation (right panel) in the same simulation. Since the 
incoming and outgoing walls are relativistic and there is little scalar radiation expelled in the collision, Eq. [8] should 
be a good estimate of the mass parameter in the region immediately to the future of the collision. Using the values 



of H in the true and false vacuum, we obtain Mm = 0.02 as an estimate. In the left panel of Fig. 14 we overplot 
(red-solid and blue dashed lines) the predicted behavior of a as a function of z; the agreement between the numerical 
solution and the thin- wall solution is excellent. A more detailed thin- wall construction takes into account the scalar 
radiation emitted by the collision. In this case, there are two mass parameters: one describing the metric in the false 
vacuum pocket and one describing the metric in the future of an outgoing null shell of scalar radiation. Assuming the 
metric in the future of the radiation shell is of the form Eq. [2] we can find the best-fit to the mass parameter for the 
simulated collision. Doing so, we find ~ 0.003; the best-fit is overplotted (blue dashed line) on the curve obtained 



from the simulation in the right panel of Fig. 14 The energy density in the outgoing shell of radiation is related to 
the mass parameter by 

(59) 



47rz' 



We can therefore estimate the initial energy density in scalar radiation released during the collision in this example to 
be tTr(zc) — 10~^i/|,. Repeating this exercise for bubbles with several different initial separations, we find that Eq. [s] 
predicts the mass parameter in the post-collision false vacuum pocket very well. The mass parameter associated with 
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FIG. 14: The metric function a evaluated at x = 0.375 (left) and x = 0.280 (right) for the collision between two identical 
bubbles formed in the potential VI and initially separated by Xc = 0.75. The plot on the left cuts through the region that 
underwent a transition back to the false vacuum. We overplot the expectation for a from the thin-wall approximation (blue 
dashed line) as well as the metric in empty hyperbolic de Sitter outside (solid red line) and inside (dot-dashed purple line) the 
bubble. The agreement between the thin-wall solution to the immediate future of the collision (at z ~ 0.4) and the numerics is 
excellent. The plot on the right cuts through the region to the future of the outgoing scalar radiation, but outside the pocket 
of false vacuum. We find the best-fit for the mass parameter, assuming that a = Aj (see Eq. bl in this region, obtaining 
Mr ~ 3.0 X 10"^. This function is over-plotted on the simulation as the blue dashed line, and in aodition we show the metric 
function for empty HdS inside (dot-dashed purple line) and outside (solid red line) the bubble. 



the region to the future of the outgomg radiation is roughly 1/10 of the mass parameter associated with the false 
vacuum pocket. Therefore, we conclude that most of the energy released in the collision goes into the formation of 
the false vacuum pocket as opposed to the outgoing scalar radiation. 

Finally, we discuss bubble collisions for the potentials V3 and V4. As described above, when two highly relativistic 
bubbles collide, the field configurations in the vicinity of the collision simply superpose. If two identical vacuum 
bubbles collide, and there happens to be a third vacuum located at a distance in field space equal to twice the 
separation between the true and false vacua, then the collision can result in a lasting region of the third vacuum. 



This phenomenon is known as a classical transition, and was introduced in Refs. |17| and |18j . In Fig. 15 we plot 
two constant-z slices through the collision of two bubbles formed from the potential V3. Just after the collision, the 
incoming field profiles superpose. In this potential, there is another vacuum of lower energy density at this position 
in field space. Therefore, a classical transition results, as shown in the contour plot in the left panel of Fig. [16 These 
solutions are in good qualitative agreement with previous simulations which neglected gravitational effects [T^l HTJ [TB] . 

In Ref. [22], the effect of gravity on classical transitions was discussed using the junction condition formalism. In 
the absence of gravitational effects, classical transitions can only lead to lasting regions of new vacua that have lower 
energy density than the interior of the colliding bubbles. However, the gravitationally repulsive nature of domain 
walls, and the background expansion of space in hyperbolic de Sitter, can allow for the formation of lasting regions 
that have higher energy density than the colliding bubble interiors, but no higher than the false vacuum from which the 
bubbles were formed. The easiest way to ensure a lasting region of new vacuum is to enclose it with domain walls that 
have enough energy density to be gravitationally repulsive (see the discussion surrounding Eq. |4|: \H'j^ — Hl\ < kf. In 
the right panel of Fig. |16[ we show a simulated collision in the potential V4 which produces a lasting region of higher 
vacuum through the production of repulsive walls. It is also possible to create a lasting region by creating marginally 
repulsive domain walls (which have a tension close to the bound for repulsive walls) with sufficient kinetic energy to 
grow a region of new vacuum larger than the size of the interior causal horizon. In practice, for the examples we 
have studied, the required initial separation of the bubbles is larger than what we can access reliably with our current 
numerical implementation. 

In summary, we find that most of the energy released in the collision of identical vacuum bubbles goes into localized 
field configurations at the position of the collision. Very little energy produced in the collision is radiated away, at 
least initially. In addition, the thin-wall solutions are in good agreement, both qualitative and quantitative, with the 
simulated bubble collisions. We have confirmed the existence of classical transitions in the presence of gravity, and 
found that lasting regions of higher energy density can indeed be produced through gravitational effects as predicted 
by Ref. [12] . 
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FIG. 15: The field on two constant-time slices (left and right) just before and just after the collision of two identical bubbles 
formed in the potential V3. The field linearly superposes in the vicinity of the collision. This value of the field is in the vicinity 
of a third vacuum, inducing a classical transition. 





FIG. 16: Contour plots of the field for the collision of two identical vacuum bubbles formed in the potential V3 (left) and V4 
(right). In the right panel, the boundaries of the simulation lie at the center of each bubble. In both cases, the collision causes 
a classical transition to a region containing another vacuum. On the left, the new vacuum has lower energy than either bubble 
interior, and therefore the region that underwent the classical transition expands. On the right, the new vacuum has higher 
energy than either bubble interior, but because the domain walls surrounding the new vacuum are repulsive, the region that 
underwent the classical transition survives. 



Non-identical vacuum bubbles 



We now move on to the case where the colliding vacuum bubbles are not identical. Since the interiors of the colliding 
bubbles are not identical, a domain wall must form in the future of the collision to separate the two bubbles. In the 
junction condition formalism, the post-collision domain wall is characterized only by its tension. To conserve energy 
and momentum at the collision, a shell of outgoing scalar radiation is assumed to be produced as well. The late-time 
behavior of the post-collision domain wall in the junction condition formalism is fully determined by its tension, and 
the vacuum energies on either side by Eq. |4] The wall can either be normal (accelerating towards the side with a 
lower vacuum energy) or repulsive (accelerating away from both sides). In this section, we confirm these late-time 
properties of the solution with our numerics, and show that taking into account the full dynamics of the collision can 
lead to far more varied behavior in the vicinity of the collision. 

As a representative example, we show a contour plot of Lp in the left panel of Fig. [17] for the collision of the two 
dilTerent types of bubbles allowed by the potential VI. Focusing for the moment on the behavior of the field in the 
immediate aftermath of the collision, note that the post-collision domain wall has a great deal of structure. A portion 
of the energy released in the collision goes into a breathing mode of the wall: a pocket of the false vacuum is created 
during the collision, which then expands and contracts. In fact, the excitation of this mode can easily be predicted 
using the free-passage approximation discussed earlier. The potential barriers Tl and T2 for the potential VI are 
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FIG. 17: On the left, we show a contour plot of (p{x,z) for the collision between the two different bubbles allowed by the 
potential VI. On the right, the potential barrier in the region Tl has been compressed by a factor of two relative to the 
potential VI. In the free-passage approximation, the structure of the potential determines the position of the field in the 
immediate aftermath of the collision. For the collision on the left, the potential barriers are nearly symmetric, and the sum 
of the wall profiles is close to ip = 0. This leads to the production of an initially expanding pocket of false vacuum inside the 
post-collision domain wall. The collision on the right was generated from a potential where the width of Tl is roughly half that 
of T2. The sum of the wall profiles is at positive i/p near the vacuum of the Observation bubble. Therefore, the post-collision 
domain wall is formed moving away from the Observation bubble interior, and no internal modes of the post-collision domain 
wall are excited. In both cases at late times the post-collision domain wall accelerates from the Observation bubble interior, as 
predicted by the junction condition formalism in Eq. [4] 

nearly symmetric. This means that when the incoming field profiles making up the bubble walls superpose, the field 
is placed very near the false vacuum at = immediately after the collision. Thus, we excite the breathing mode 
of the post-collision wall. Such modes increase the energy density and tension carried by the wall, and change the 
dynamics, causing the wall to accelerate more slowly out of the Observation bubble. 

The free-passage approximation predicts that fewer internal modes of the wall are produced when the potential 
barriers making up the colliding bubbles are asymmetric. To test this, we generate a potential identical to VI, but 
where the width of the barrier Tl is halved. Superposing the bubble wall profiles, the field should be near the vacuum 
inside the Observation bubble immediately after the collision. A contour plot of ip for the collision in this case is 
shown in the right panel of Fig. |17[ Here, the wall forms already moving away from the Observation bubble interior 
with no breathing modes, as expected. If the width of the barrier Tl is increased, then superposing the bubble wall 
profiles, the wall is predicted to form moving into the Observation bubble with no breathing modes. This is again 
observed in our simulations. 

After a short period of time, the breathing modes of the wall damp, and for the examples shown in Fig. |17[ the wall 
accelerates out of the Observation bubble. How does this compare to our expectation from the junction condition 
solutions? From the potential, the difference in Hubble parameters is Hj^ — iJ| ~ 0.13, and we can estimate the 
tension of the post-collision domain wall for the potential VI as the sum of the pre-coUision domain walls, yielding 
~ 0.001. Since the difference in Hubble parameters is much greater than the tension (making this domain wall 
of the normal type), we expect that at late times the post-collision domain wall accelerates in the direction of the 
bubble with lower vacuum energy: in this case, the Collision bubble. This is indeed what we observe. Again from 
the junction condition solutions, we expect that the potential V2 gives rise to a repulsive post-collision domain wall. 
This is because, while the tension of the post-collision domain wall is approximately the same as for VI, the difference 
in Hubble parameters is much smaller: iJ^ — iJ| ~ 0.0006. We plot this example in Fig. [g] Calculating the proper 
distance from the center of each bubble to the wall, we find that it is indeed accelerating away from both bubble 
interiors. 

In summary, we have confirmed that the late-time behavior of the domain wall produced in the collision between two 
different vacuum bubbles is well described by the junction condition formalism, but also discovered a rich dynamics in 
the immediate aftermath of the collision. For the potentials we have studied, the relative widths of the barriers Tl and 
T2 determine the outcome of the collision, as predicted by the free passage approximation. For symmetric barriers, 
the collision excites a breathing mode of the post-collision domain wall. For asymmetric barriers, the post-collision 
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domain wall forms already moving in the direction of the thicker barrier, with no breathing modes. 



C. Colliding bubbles with an interior cosmology 



We now turn to the case where there is an inflationary potential inside the Observation and Collision bubbles. In 
this case, it is not possible to apply the Israel junction condition formalism to predict the outcome of bubble collisions: 
the resulting spacetime depends on the full non-linear dynamics. This section is also most relevant for determining 
the possible observational signatures of bubble collisions. As we discussed previously, our current implementation is 
not well-suited to run for the timescales necessary to extract the signal accurately. However, we can gain much insight 
into the dynamics of the collision, and the expected amplitude of any possible signals, by studying the outcome of 
the collision over a limited number of e-folds of the inner-bubble cosmology following the collision. 

As a flrst application, we collide two identical bubbles formed from the "large- field" potential LI. A contour plot of 



the field is shown in Fig. 18 In the case of large-field models, it is necessary to design potentials with a hierarchy of 
scales between the width of the barrier T2 and the width of the inflationary segment C2. This is evident in Fig. [5] and 
is necessary because potentials with a barrier width comparable to Mpi typically do not support CDL bubbles, while 
inflation occurs as the field rolls over a distance greater than Mpi (see Ref. p[2 fo'" further discussion of this point). 
This is quite relevant for determining the outcome of the collision, because in the free passage approximation, the field 
is moved a distance comparable to the barrier width down the inflationary potential. This distance is predicted to be 
relatively small compared to the total field excursion in large-field models. In Fig. [19] we plot two constant-z slices 
through the simulation just before and just after the collision. The red dashed line represents the sum of the incoming 
wall profiles, and it can be seen that initially the field profiles do indeed simply superpose. After the collision, inflation 
proceeds, although from a different initial value of the fleld as the rest of the bubble interior. We confirm that infiation 
is still occurring in this region by checking that the metric functions are still approximately HdS in the future of the 
collision. Examining the colorbar in Fig. 18 and Fig. [s] the field has not progressed very far down the potential. The 
conclusion is that there are fewer e-folds of inflation in the future of the collision than in the undisturbed portions of 
the bubble. We expect that this is a generic prediction for the outcome of the collision between identical bubbles. 

We have also examined collisions between bubbles with different interior cosmologies generated from the potentials 
LI and L2. These potentials differ only in the choice of Tl and CI: for LI, the potential slopes down from the barrier 
Tl to a vacuum, while for L2, there is a potential minimum immediately on the other side of the barrier Tl. Using the 
potential LI, we can simulate the collision between two bubbles with an interior cosmology, while using the potential 
L2, we can simulate the collision between a vacuum bubble and a bubble containing an interior cosmology. 



In Fig. 20 we show a simulated collision between the two bubbles formed from the potential L2. Since in this case 
the energy is lower everywhere inside the Observation bubble, we expect the post-collision domain wall to accelerate 
out of the Observation bubble at large z; this is conflrmed in the simulation. Since the potential barriers are nearly 
symmetric in this case, the free-passage approximation predicts that internal modes of the post-collision domain 
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FIG. 19: The field on two constant-z slices just before and just after the collision in the simulation shown in Fig. |18| The 
red-dashed hue is the sum of the distance in field space between the false vacuum and the position the endpoint of the CDL 
instanton. In agreement with the free passage approximation, the field profiles in the immediate future of the collision simply 
superpose. 




FIG. 20: A contour plot of f for the collision of two different bubbles formed from the potential L2. In the immediate future 
of the collision event, the inflaton inside the Observation bubble is perturbed, and a pocket of the false vacuum is formed. 
This pocket of false vacuum is a breathing mode of the post-collision domain wall, which once formed, accelerates out of 
the Observation bubble. The breathing mode further perturbs the inflaton inside the Observation bubble. The amplitude of 
perturbations to the inflaton is small compared to the total field excursion during inflation in this model. 

wall are excited; this is also observed in the simulation. In the future of the collision, the field is pushed down the 
inflationary plateau, getting an extra kick each time the post-collision domain wall undergoes an oscillation. The 
amplitude of these perturbations is somewhat smaller than it was for the collision between identical bubbles. This is 
to be expected, as the energy of the collision is transferred to the creation of the post-collision domain wall and the 
excitation of breathing modes. Increasing the center of mass energy of the collision by increasing the initial bubble 
separation mainly excites larger amplitude breathing modes, as opposed to larger perturbations of the inflaton inside 
the Observation bubble. This is consistent with the results of Sec. |VI B[ where it was found that for vacuum bubbles 
very little energy is released in scalar radiation after the collision over a wide range of kinematics. 

As in the case of colliding vacuum bubbles, we observe that fewer internal modes of the post-collision domain wall 
are excited when the potential barriers making up the two different bubbles are asymmetric. Again, we adjust the 
width of the Tl barrier while keeping the other properties of the potential L2 constant: in one case stretching it by 
a factor of two and in another case compressing it by a factor of two. A contour plot of the field is shown in Fig. [2l] 
for the thick barrier Tl on the left, and the thin barrier Tl on the right. In the case where Tl is twice as thick 
as T2, the inflaton inside the Observation bubble is hardly perturbed at all! This corroborates the existence of the 
"mild-collisions" first discussed in Ref. [TT]. In the case where Tl is half as thick as T2, the perturbation is somewhat 
smaller than for symmetric potential barriers. 
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FIG. 21: On the left, we show a contour plot of Lp for the collision of bubbles in a potential with a barrier Tl twice as large as 
that of the potential L2. Initially, the post-collision domain wall moves into the Observation bubble, but then turns around. 
There is very little disturbance to the inflaton inside the Observation bubble. On the right, we show a contour plot of ip for 
the collision of bubbles in a potential with a barrier T\ half as large as that of the potential L2. The post-collision domain 
wall forms accelerating out of the Observation bubble, and no breathing modes are excited. These two examples are consistent 
with the predictions of the free passage approximation. In both panels, the boundaries of the simulation lie at the center of 
each bubble. 

We have studied a number of different initial separations for the bubbles, finding that the field excursion in the 
aftermath of the collision is bounded by the sum of the field profiles; this confirms the free-passage approximation is 
valid. What changes with increasing bubble separation is the size of the spacetime region inside which the profiles 
superpose. Thus, once the bubbles are sufficiently relativistic, increasing the energy in the collision only goes into 
creating larger regions inside which the field is displaced a constant distance. This suggests that the amplitude of the 
observable signature of bubble collisions is more sensitive to the structure of the potential than to the kinematics of 
the collision. 

We now add a cosmology inside the Collision bubble, simulating collisions in the potential LI. This can drastically 
affect the outcome of a bubble collision, as shown in Fig. [22j In this example, there is a relatively steep potential 
slope inside the Collision bubble. Even though the energy inside the Collision bubble is initially higher than inside the 
Observation bubble, the collision causes the field to quickly roll to a much lower energy than the inflationary plateau 
inside the Observation bubble. Thus, the post-collision domain wall begins to eat into the Observation bubble. After 
inflation, since the ultimate vacuum in the Observation bubble is much lower than the vacuum inside the Collision 
bubble, we expect that the wall will retreat and move into the Coflision bubble. Our current implementation cannot 
reliably track the evolution for such long timescales, and so we cannot conflrm this numerically. 

In Ref. [12], it was suggested that inflationary potentials of the "small- field" type would be more susceptible to 
disruption by collisions than potentials of the large-field type studied above. However, since this work did not include 
gravitational effects, this was not conclusively established. The free-passage approximation discussed above already 
provides some idea of which models of inflation are susceptible to disruption by a collision. Focusing for the moment 
on a collision between identical bubbles, the field in the immediate aftermath of the collision is displaced a distance 
equal to the width of the barrier making up the bubble. Models of infiation which have an inflationary plateau smaller 
than or comparable to this distance are susceptible to complete disruption by a collision. This implies that there 
should be a hierarchy in scales between the width of the barrier T2 and the inflationary potential C2 for inflation to 
survive bubble collisions. 

The two types of small-field inflation we use in our numerics are shown in Fig. [5] For the potential SI, inflation 
occurs near an inflection point (see e.g. Refs. [S71[5S] for a discussion of similar models). For the potential S2, inflation 
occurs near a potential maximum (see e.g. Ref. |59|). The inflaton has a much smaller range in SI than S2, although 
in both cases, the held excursion during inflation is sub-Planckian. 

In Fig. [23] we show contour plots for the collision of bubbles in the potential SI. On the top left is the case where 
the colliding bubbles are identical. In agreement with the non-gravitational simulations of Ref. [T^j, the collision 
completely ends inflation inside its future light cone. This is expected since the width of the barrier Tl is much 
wider than the region of the potential C2 where the slow-roll parameters are satisfled. In the top right panel, we 
adjust the barrier width T\ to be one fourth the width of the barrier T2. Again, inflation is completely disrupted by 
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FIG. 22: A contour plot of ^) for the collision of bubbles in the potential LI. Because the slope CI away from the potential 
barrier is much steeper inside the Collision bubble than inside the Observation bubble, the post-collision domain wall moves 
into the Observation bubble during its inflationary phase. 

the collision. However, adjusting the barrier width Tl to be four times the width of the barrier T2, the inflaton is 
hardly disturbed, as shown in the bottom panel. We found a similar effect for the large-field models discussed above. 
Therefore, the sensitivity of inflection point models of inflation to bubble collisions depends quite sensitively on the 
structure of the potential barriers making up the colliding bubbles. 

We have also simulated collisions in the potential S2. In this case, the barrier width is somewhat smaller than 
the field excursion during inflation, and we expect that collisions will not completely disrupt inflation inside the 
Observation bubble. This is indeed what we observe in our simulations. The qualitative picture is nearly identical to 
the large-field cases discussed above. However, the duration of inflation inside the region affected by the collision will 
be somewhat less than it was in the large-field examples. 

Let us summarize the results of this section. For the collision between identical bubbles, the inflaton is displaced 
down the inflationary potential a distance comparable to the barrier width, in agreement with the free passage 
approximation. Inflation has a shorter duration to the future of the collision than in the undisturbed portions of 
the bubble. In large-field models, this displacement is rather small compared to the total field excursion. In small- 
field models, the displacement can be so large that inflation is completely disrupted. Increasing the initial bubble 
separation appears to have a relatively small effect on the number of e-folds in the future of the collision. For non- 
identical bubbles, the structure of the potential barriers is the dominant factor in determining the outcome of a 
collision. This is again in agreement with the free passage approximation. When the barrier is symmetric, breathing 
modes of the wall are excited. When the potential barrier Tl is thinner than the barrier r2, the post-collision domain 
wall is formed moving out of the Observation bubble without breathing modes. When the potential barrier Tl is 
thicker than the barrier T2, the inflaton inside the Observation bubble is minimally disturbed. Even the inflection 
point model of inflation (potential SI) can survive a collision in this case. If there is a cosmology in both colliding 
bubbles, the late time behavior of the post-collision domain wall depends on the details of the field evolution. For 
example, in the potential LI the energy density in the Collision bubble quickly falls below that of the Observation 
bubble (which is undergoing slow- roll), causing the post-collision domain wall to have multiple turning points. For 
the post-collision domain wall to monotonically progress out of the Observation bubble, the Collision bubble must be 
undergoing slow-roll at an energy higher than that of the Observation bubble. Finally, we note that in each of the 
models we have studied, there are always fewer e-folds of inflation to the future of the collision than in the undisturbed 
portions of the bubble. 

VII. CONCLUSIONS 

This work establishes the groundwork for rigorously determining the connection between an underlying potential 
landscape giving rise to eternal inflation and the cosniological signal of bubble collisions. Such a goal requires the 
ability to study the coupled Einstein and scalar field equations without approximations. To this end we have developed 
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FIG. 23: Contour plots of ip for the collision of bubbles in the potential SI. On the top left, we simulate the collision between 
two identical bubbles. In the future light cone of the collision, inflation is completely disrupted (as can be seen by comparing 
the color bar to the field values on the potential in Fig. [5|. On the top right, we simulate the collision of two different bubbles 
in the case where the width of the barrier Tl is adjusted to be a quarter of the width of the barrier T2. Again, inflation is 
completely disrupted by the collision. On the bottom, we simulate the collision of two different bubbles in the case where the 
width of the barrier Tl is adjusted to be four times the width of the barrier T2. In this case, the inflaton inside the Observation 
bubble is hardly perturbed, and inflation continues to the future of the collision. In all cases, the boundaries of the simulation 
lie at the center of each bubble. 



a robust code that implements the full system of relativistic field equations, thoroughly tested it on a number of cases 
and applied it to determine the outcome of collisions between bubbles. 

Our results highlight the utility of two complementary analytic tools for predicting the outcome of bubble collisions. 
For the collision of vacuum bubbles, we have confirmed that the thin-wall solutions obtained from the Israel junction 
condition formalism approximate the numerical solutions to good accuracy. In particular, the junction condition 
solutions give a good quantitative prediction for the metric functions to the future of the collision, as well as a 
good qualitative prediction for the late-time behavior of the post-collision domain wall. The limitation is that the 
junction condition formalism cannot capture the dynamics of the collision event. This is particularly problematic 
when considering bubbles with an interior cosmology. The dynamics are, however, captured to a certain extent by the 
free passage approximation [T5] which states that, in the immediate future of the collision, the field profiles making 
up the two bubbles simply superpose. This illustrates that the structure of the potential can be more important for 
determining the outcome of a collision than the kinematics. The free passage approximation breaks down shortly 
after the collision, at which point the full field dynamics must be considered. A full numerical solution is necessary to 
bridge the gap between the realm of validity of the free passage approximation and the junction condition formalism. 

Some of the main conclusions we have reached based on our numerical solutions are as follows: 

• The energy released in the collision of identical vacuum bubbles goes mostly into the formation localized field 
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configurations such as oscillons and pockets of the false vacuum. 

• Classical transitions occur in the presence of gravity, and lasting regions of higher energy density can be created 
through gravitational effects. 

• The structure of the potential is the dominant factor in determining the immediate outcome of a collision. 
For the collision between identical bubbles, the field is initially displaced a distance equal to the barrier width 
towards its minimum. For the collision of non-identical bubbles, when the potential barriers are symmetric, the 
collision excites a breathing mode of the post-collision domain wall. For asymmetric barriers, the post-collision 
domain wall forms with some initial velocity in the direction of the thicker barrier, with no breathing modes. 

• We have demonstrated conclusively that slow-roll infiation can occur to the future of a collision. Large field 
models of infiation are more robust to collisions than small-field models, which are completely disrupted unless 
the potential barriers are highly asymmetric (with the barrier Tl much thicker than the barrier T2). 

• The direction in which the post-collision wall accelerates is determined not only by the height of the potential 
at the exit from the barrier, but also by the slope of the potential as the field tends towards the minimum. 
Therefore, if there is an inflationary potential in the Observation bubble, for the Collision bubble not to intrude, 
the interior of the latter is required to be either a vacuum (no cosmology) or an inflationary plateau that is 
either higher or broader than the one in the Observation bubble. 

• The amplitude of the perturbation to the inflaton is set mostly by the barrier widths in the potential, with only 
a mild dependence on the kinematics. In addition, for the examples we have studied, there are always fewer 
e-folds to the future of the collision than in the undisturbed portions of the bubble. 

Based on these observations, we can determine which potentials are most likely to give rise to cosmologies including 
bubble collisions that are compatible with our current observable universe, yet leave detectable signatures. Potentials 
resembling L2 in Fig. [5] are good candidates. Since this potential includes a large- field model of infiation, the collision 
does not disrupt cosmological evolution inside the bubble. Adjusting the relative width of the two barriers changes 
the amplitude of the perturbation to the inflaton. The intrinsic amplitude of the perturbation along with the total 
duration of infiation inside the Observation bubble determines the amplitude of the signal in the CMB or other 
cosmological observables. 

While the existing code is able to accurately obtain the dynamics of the collision, it cannot yet be used for 
studying the ensuing cosmology: the uniform grid employed and coordinate conditions adopted are unable to track 
the increasingly narrow structures that develop around walls in a computationally efficient manner. To address this 
issue, we will extend our code to adopt adaptive mesh refinement (see e.g. Refs. [iQl [60] for a first example of 
its application in numerical relativity). We will present the derivation of the observational signatures from these 
simulations in a forthcoming publication. 
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Appendix A: General equations of motion 



In this appendix, we derive the general form of the equations of motion assuming the line element Eq. [9j In this 
appendix, we work in units where Mpi = 1. The results of this appendix are used in Sec. Ill to derive the gauge- fixed 
equations of motion used in the simulations. 
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a. Stress energy tensor and constraints 



We begin by finding the energy momentum tensor and field equations for the scalar field. The energy momentum 
tensor for the scalar field is given by: 



The non-trivial components of this tensor are, 
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The density, current and spatial stress associated with T'^" determined by a normal observer are given by, 
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The Hamiltonian and momentum constraint equations are given by: 
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which for our problem of interest reduce to 
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(only the x component of the momentum constraint is non-trivial) . 



6. Evolution equations 



The evolution equation for the scalar field is given by 
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This yields the evolution equation for 11, 
n : 
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The evolution equation for $ is obtained by taking the spatial derivative of ip: 

$ = 9, (/3ci> + ^n) , 

where, in terms of our variables. 
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completely determine the equations of motion for the scalar field Lp. 

The evolution equations for the three-metric are obtained from the definition of the extrinsic curvature, 
and X — X components yield: 

a = -aaK'^x + /3a' + a^' , 
6 = -ahK^^ +I3b' 

Einstein equations determine the evolution equations for the extrinsic curvature, 
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Appendix B: Recovering hyperbolic dS 
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The system of equations Eq. 29 30 31 and 32 should recover the hyperbolic foliation of dS space in the case where 
<f> = n = and V — const. In this appendix, we work in units where Mp\ ~ 1. The metric in hyperbolic dS is given 
by: 



ds^ = -(1 + H^z^ydz^ + (1 + H^z^)dx^ + z^dH^ . 



The evolution equation for 6, Eq. (15), defines 
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This expression can be replaced in the Hamiltonian constraint 

2 



167rF(^) = + ^K^'^K'^^ + 2ifx^2 ^ 



(Bl) 



(B2) 



(B3) 



30 



to determine an equation for ii'^^ 



K\ = - . (B4) 



This can now be exploited by noting the momentum constraint, 

d^K\ = 0, (B5) 
implies that a — a{z). One can then find a{z) by solving the evolution equation for K^^: 

i.x,= 3^11±^^:M^. (B6) 



Exploiting again Eq. (B2|, one obtains an equation for a: 
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The general solution of this equation is: 
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a=(^l + — z^-— j , (B8) 

where we have set an integration constant to M. For the case M = 0, and identifying = 87ry((^)/3, this is the 
metric function we expected. 

The evolution equation for a yields: 

a = -otaK''^ = ~ + ao?z\-KV{if) . (B9) 

Substituting now with the expression obtained for a, 

a= y^', a, (BIO) 

whose solution is given by: 
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which is again the expected solution. 



Appendix C: Field equations and characteristic structure in hyperbolic de Sitter 

The equation of motions for a scalar field configuration possessing hyperbolic symmetry in a background HdS 
spacetime is given, in terms of the first order variables {ip, 11, $} by: 

ip = {1 + H^z^)-^U (CI) 

n = {i + H^z^}-^^' --n-d^v{ip) (C2) 
$ = (1 + i^2z2)-l^' . (C3) 

having identified — 8ttV/3. The eigenvalues of the principal part of these equations determine the characteristic 
speeds of the system, which are given by 

c = ±{l + H^z^)-^;0. (C4) 

Thus, the physical (non-zero) characteristic speeds range from unity at small z to approaching zero as z becomes 
large compared to H. 
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V[x], F„[x] 




-1.5 -1.0 -0.5 0.5 1.0 1.5 

FIG. 24: V (black) and vo (red dashed) with a = 0.1 and C = 2 



Appendix D: Numerical algorithm for finding CDL instantons 



In this appendix, we outline a numerical method for determining the initial data for single bubbles. The initial data 
for a single bubble spacetime is given by the analytic continuation of the Coleman-de Luccia (CDL) instanton [2]. 
The CDL instanton is a solution to the Euclidean Einstein and field equations that describes the formation of a 
bubble. First, we outline a convenient set of dimensionless variables and the associated equations of motion. Then, 
we describe our numerical algorithm for finding the CDL instanton. Finally, as an aside, we apply our numerics to 
test the validity of the thin-wall approximation. 

1. Variables and equations of motion 

Consider the class of potentials defined by 

V{^) = /i; i^/M) , (Dl) 

We define the following dimensionless variables: 



M'^^li"'*"!/'*^ " 3M? 



PI 



Substituting these into the Euclidean Einstein (Eq. [35|) and field (Eq. 37 ) equations, we obtain: 



3f 

x+ —x-v' ^0, (D3) 
r 



(±2 + v) , (D4) 



where the primes and dots, respectively, refer to x~ and s— derivatives. To find the CDL instanton, we must solve a 
double-boundary value problem with initial and final conditions given by: 

X{s = 0) ~ XT, x{s = 0) = 0, r(s = 0) = 0, x{s = Smax) ^ Xp, X{s = Smax) = 0, r{s = Smax) = . (D5) 

Such solutions interpolate between the basin of attraction of the true and false vacuum as the scale factor r goes 
between its two zeroes at s = and s — Smax- An example is shown in Fig. [2j 
In this appendix, we consider potentials of the form 

X X „ ax^ 
where C and a are constants. An example is shown in Fig. |24[ 
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Since the instanton solution interpolates between the false and true vacua in the classically forbidden region of the 
potential, we can use this solution to obtain the probability per unit four-volume that a bubble forms from the false 
vacuum in the WKB approximation. This is given by 0: 

r = Ae-^, (D7) 

where 

B = Si-Sbg- (D8) 

The first term is the Euclidean action, evaluated for the instanton solution. Substituting the equations of motion into 
the Euclidean action for a canonically normalized scalar field minimally coupled to gravity, the instanton action is 
given by 



S, = -47r2 



Af4 

M / Js=0 



ds(-r^v+'^^ . (D9) 



where r(s) is the instanton solution. The second term in Eq. D8 is known as the background subtraction, and is 
simply the Euclidean action of the false vacuum de Sitter space: 

Sbg^~^. (DIO) 



2. Numerical strategy 



To solve for the CDL instanton, we utilize the over/under-shoot argument of Coleman. We have formulated a 
numerical implementation of the algorithm described in Ref. [41] . which involves the following steps: 

1. Start the field at xq ^ xt — 10~^^. We start the numerical evolution aX s = ds to avoid the singular nature 
of the friction term at r = 0. The initial conditions for the evolution ai s = ds can be solved for by Taylor 
expanding the solution in s near s — Q. The first non-trivial dependence on the potential arises at 2nd and 3rd 
order respectively for x and r. Expanding to 3rd order, and solving for the coefficients in the Taylor expansion 
(the first and third-order terms in x are zero; the 2nd order term in r is zero), we obtain: 



x{ds) = xa-\- 



dxv{xo)ds'^ 



x(ds) 



dxv{xo)ds 



ids) 



^v{xo)ds^ 
6 



(Dll) 



We then evolve the field and test for an overshoot. 
2. If an overshoot is found on the original test: 

(a) Define xl = a^max and xr = xt, then bisect and set xq equal to the midpoint: xq = xl + {xr — xl)/'^- If 
there is an overshoot, xq is too close to xr, so redefine xr = Xq and test again. If there is an undershoot, 
Xq is too close to xl, so redefine xl = xq and test again. Repeat these tests until {xr — Xl)/'^ is at a 
pre-specified tolerance (we use 10~^*, which is of order machine precision). In doing so, we iterate to the 
solution of the instanton end-point in the vicinity of the true vacuum. 

(b) The evolution is terminated either when x — xp or x changes sign. Either way, the scale factor has not 
yet reached its second zero. Therefore, we match onto a solution where: 



x{s 



end } 



r{s) 



1 



i^v{xf) 



z^Jv{xf){s - s*) 



with 



Send — arcsm 



\/v{xF)r{scr,d) /{e^/v{xF)), r>0 



■Send — TT -t- arcsm e 



\/v{xF)r{scnd) /{t\/v{xF)), r < 



(D12) 

(D13) 
(D14) 



where we have been careful to take the correct root of the arcsin function. This solution is used until r = 0. 
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3. If an overshoot is not found on the original test, the field must be too close to the true vacuum for us to resolve 
numerically. 

(a) We use the solution 



x = xt, r{s) = — sin e\/ v{xt)s (D15) 

ey/v{xT) L J 

for a time Asi, substitute in the new initial conditions for x and r at that time, and test for overshoot. 
If no overshoot is found, use the solution for another period Asi, and re-test for overshoot. Repeat until 
overshoot is found. 

(b) Now, find the instanton end-point by bisection as above. 

(c) As above, find the portion of the solution near xp- 

4. Calculate the instanton action and export the initial data. 
The output of the instanton calculation is the post-tunnelling field configuration: 

^4^. (D16) 
We then convert the solution into units appropriate to the simulation: 



^ Mpi M V Stt j M ■ ^^^'^^ 
Additionally, we have 

i = ey/v{(pF) s (D18) 

V=^4^,. (D19) 



3. Checking the thin-wall approximation 

Potentials where the field loiters in a close neighborhood of the true vacuum for a time much longer than the 
Euclidean transition time from the true to the false vacuum are known to satisfy the "thin- wall" approximation [5]. 
In this approximation, the bubble wall is assumed to be a rigid membrane with tension a. The instanton can be 
constructed by pasting a portion of a true vacuum dS 4-sphere to a portion of a false vacuum dS 4-sphere. The 
tension accounts for the discontinuity in the derivative of the metric across this junction. We can check our numerical 
solutions against the thin- wall approximation. In this limit, equations for the critical radius of the bubble and 
tunnelling exponent were first derived in Ref. [61], and are given by 



9a2 



crit 



{v{xf) - v{xt)) 



3ea{v{xFy-v{xT))^ 
^ {v{xf)+v{xt)) 



3e(T 



2{v{xf) -v{xT)yr^ 



(D20) 



and 



where 



B = 



2{v{xf) ~ vixrW 



r{x,y) 



Sea 



v{xf) + v{xt) 
2{v{xf) — v(xT)y^'^ J 'v{xf) — v{xt) 



(1 + xy) - + 2xy + x'^ 



'■{y^ - + 2xy + x"^ 



(D21) 



(D22) 



The vacuum energy is determined by the potential, and therefore to evaluate the thin-wall formulas above we need 
only calculate the tension, defined by: 



(7 = 2/ ds{vo[(j){s)] - Vo[(I)f]) : 



(D23) 
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B 



SxlO' 

( 

2x10' 
IxlO' 
SxlO"* 

2x10" 
IxlO" 



0.002 



0.005 0.010 0.020 



0.050 0.100 



FIG. 25: The tunnelling exponent as a function of e computed for a few examples numerically (dots) and from Eg. |D2l] using 
the upper and lower estimates of the tension (red-dashed lines). 



Si, Sbg 

10" 
10" 

10" 

10^ 



0.001 0.002 0.005 010 0.020 0.050 100^ 

FIG. 26: The numerically computed instanton action (dots) plotted against the background subtraction (red-dashed line). 
Both are growing like e~* as e — >■ 0. 



where vq is a potential without the symmetry-breaking terms. For potentials of the form Eq. D6 vq is just 

„2 ^4 



X 



c. 



(D24) 



This is shown in Fig. [24j We obtain slightly different estimates for the tension by changing the limits of integration 
(the numerical values we obtain are ctiow = 0.94 and (Thigh = 1-02). The numerical value of the tunnelling exponent is 
quite sensitive to any errors, since it depends on the fourth power of the tension. 



In Fig. 25 , we compare the tunnelling exponent obtained from Eq. D21 to that obtained by integrating the numer- 
ically determined instanton solution (using Eq |D9| and Eq. DIO) as a function of e. The data points follow the curve 
well for large values o f e, but then fall off for smaller e due to truncation error. In Fig. |26[ the numerically computed 
instanton action (Eq D9) and the background subtraction are plotted separately. Both grow like as e — > as 
expected, and must cancel increasingly more precisely as e goes to zero (where the instanton action should approach 
a constant). 

Our numerical solutions corroborate the validity of the thin-wall equations for the critical radius and tunnelling 
exponent. As far as the authors are aware, this is the first explicit numerical test of the validity of the thin- wall 
This 



D21 does not correctly 
all of the 



approximation. This result contradicts recent claims [621 |63] in the literature that Eq 
reproduce the instanton action Eq. |D9| computed using full solutions to the equations of motion. For 
examples we have studied, the thin-wall formulas reproduce the numerically computed tunnelling exponents to high 
accuracy. We conclude that the thin-wall approximation is indeed valid and highly accurate. 
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Appendix E: Potential Parameters 



We present here the parameters used to generate the potentials defined m Eq. 49 and shown in Figs. |4] and |5] These 
are specified in Table |l] 



type 


ffli 


a2 


M 




Vol 






Vm2 




'P32 


Ci 


u 


A 


VI 


0.5 


0.75 


3.45 X 10" 


-3 


5.0 X 


10" 


11 


N/A 


N/A 


N/A 


N/A 


2.0 


N/A 


N/A 


V2 


0.749 


0.75 


3.45 X 10" 


-3 


5.0 X 


10" 


11 


N/A 


N/A 


N/A 


N/A 


2.0 


N/A 


N/A 


V3 


-0.5 


0.50 


3.45 X 10" 


-3 


1.0 X 


10" 


12 


N/A 


N/A 


N/A 


N/A 


1.0 


N/A 


N/A 


V4 


-0.5 


1.0 X 10"® 


3.45 X 10" 


-3 


1.0 X 


10" 


12 


N/A 


N/A 


N/A 


N/A 


0.51 


N/A 


N/A 


LI 


0.5 


0.75 


3.45 X 10" 


-3 


5.0 X 


10" 


11 


-0.08 


2.75 


-0.01 


0.015 


2.0 


N/A 


N/A 


L2 


0.5 


0.75 


3.45 X 10" 


-3 


5.0 X 


10" 


11 


N/A 


2.75 


N/A 


0.015 


2.0 


N/A 


N/A 


SI 


0.5 


0.75 


3.45 X 10" 


-3 


5.0 X 


10" 


11 


-0.08 


2.75 


-0.01 


0.015 


2.0 


5.0 X 10~^ 


3.7 X 10 


S2 


0.5 


0.75 


3.45 X 10" 


-3 


5.0 X 


10" 


11 


N/A 


0.6 


N/A 


0.5 


2.0 


N/A 


N/A 



TABLE I: The parameters used to generate the potentials determined by Eq. |49| and shown in Figs. |4] and [5] 



Appendix F: Size of the simulation and resolution requirements 



Here, we discuss the maximum size of the simulation box necessary to contain two colliding bubbles. Curves of 
constant acceleration in a background HdS space take the following form when centered about 2; = 0: 



cos Xq 



;2 



(Fl) 



At large z, this asymptotes \,o x — ±.-k/2 independent of the value of xq. Now, consider the wall of the Observation 
bubble at Axq and the wall of a colliding bubble of initial radius Aaic centered about x — Xc- The two curves are 
parameterized by: 



X = arccos 



cos Ain 
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Xc + arccos 



cos Axc 



(F2) 



The maximum separation between the centers that allows for a collision occurs for Xc = tt. We can also conclude that 
the size of the box that entirely contains both bubbles is given by 



< X < Xr 



TT 

2 



(F3) 



Since the bubbles are embedded in an expanding spacetime, the resolution in terms of physical distance is decreasing 
with z. The physical distance between two points at constant z is given from the metric by: 



iJ^As = VTT^Ai 



(F4) 



At large-z, the resolution in terms of physical coordinates is decreasing roughly linearly. 

To accurately simulate bubble collisions, we must have enough resolution to describe the bubble walls at the time 
of the collision. We can estimate the required resolution as follows. If we define the wall by an inner and outer surface 
of constant field, we can make the approximation that they follow lines of constant acceleration in the background 
false vacuum de Sitter space. If these constant field surfaces are at the positions iin and Xont on the z — surface, 
they are separated by a distance 



Ax(z) 



cos Xout 



?2 



arccos 



cos Xu 



^/^ 



- (cos Xii] 

z 



cos Xout ) 



(F5) 



where the last equality holds at large z. We can therefore estimate that if TVg grid points are necessary to resolve 
the wall at small-z, we need ZcN^ grid points to resolve the wall at the time of the collision, and Zmax-^o grid points 
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to resolve the walls everywhere in the simulation. For our standard resolution of 32 x 10^ spatial grid points, there 
are typically on the order of 100 spatial grid points describing each bubble wall on the initial timeslice. Therefore, to 
retain tens of grid points in the wall at late times, runs are limited to a total duration of order z 10. 



[1] S. R. Coleman, Phys. Rev. D15, 2929 (1977). 

[2] S. R. Coleman and F. De Luccia, Phys. Rev. D21, 3305 (1980). 

[3] A. Aguirre, in Beyond the Big Bang (Springer, 2008). 

[4] A. Aguirre, M. C. Johnson, and A. Shomer, Phys. Rev. D76, 063509 (2007), arXiv:0704.3473 [hep-th]. 
[5] J. R. Gott and T. S. Statler, Phys. Lett. B136, 157 (1984). 

[6] J. Garriga, A. H. Guth, and A. Vilcnkin, Phys. Rev. D76, 123512 (2007), hep-th/0612242. 

[7] S. W. Hawking, I. G. Moss, and J. M. Stewart, Phys. Rev. D26, 2681 (1982). 

[8] Z.-C. Wu, Phys. Rev. D28, 1898 (1983). 

[9] I. G. Moss (1994), gr-qc/9405045. 
[10] B. Frcivogcl, G. T. Horowitz, and S. Sheukcr, JHEP 0705, 090 (2007), hcp-th/0703146. 
[11] A. Aguirre and M. C. Johnson, Phys. Rev. D77, 123536 (2008), 0712.3038. 
[12] A. Aguirre, M. C. Johnson, and M. Tysanner, Phys. Rev. D79, 123514 (2009). 
[13] S. Chang, M. Kleban, and T. S. Levi, JCAP 0804, 034 (2008). 
[14] S. Chang, M. Kleban, and T. S. Levi, JCAP 0904, 025 (2009), 0810.5128. 

[15] B. Czech, M. Kleban, K. Larjo, T. S. Levi, and K. Sigmdson, JCAP 1012, 023 (2010), 1006.0832. 

[16] B. Freivogcl, M. Kleban, A. Nicolis, and K. Sigurdson, JCAP 0908, 036 (2009). 

[17] R. Easther, J. T. Giblin, Jr, L. Hui, and E. A. Lim, Phys. Rev. D80, 123519 (2009), 0907.3234. 

[18] J. T. Giblin, Jr, L. Hui, E. A. Lim, and I.-S. Yang, Phys. Rev. D82, 045019 (2010), 1005.3493. 

[19] M. Kleban, T. S. Levi, and K. Sigurdson (2011), 1109.3473. 

[20] E. A. Lim and D. Simon (2011), 1103.4300. 

[21] A. Dahlen, Phys.Rev. D81, 063501 (2010), 0812.0414. 

[22] M. C. Johnson and I.-S. Yang, Phys. Rev. D82, 065023 (2010), 1005.3506. 

[23] A. Aguirre and M. C. Johnson, Rept.Prog.Phys. 74, 074901 (2011), 0908.4105. 

[24] C. L. Bennett et al. (WMAP), Astrophys. J. 583, 1 (2003), astro-ph/0301158. 

[25] E. Komatsu et al. (WMAP Collaboration), Astrophys. J. Suppl. 192, 18 (2011), 1001.4538. 

[26] S. M. Feeney M. C. Johnson, D. J. Mortlock, and H. V. Peiris, Phys.Rev. D84, 043507 (2011), 1012.3667. 

[27] S. M. Feeney M. C. Johnson, D. J. Mortlock, and H. V. Peiris, Phys.Rev.Lett. 107, 071301 (2011), 1012.1995. 

[28] L. Lchner, Classical and Quantum Gravity 18, 25 (2001), arXiv:gr-qc/0106072. 

[29] F. C. Adams, K. Freese, and L. M. Widrow, Phys. Rev. D41, 347 (1990). 

[30] J. Garriga and A. Vilenkin, Phys. Rev. D45, 3469 (1992). 

[31] J. Garriga and A. Vilenkin, Phys. Rev. D44, 1007 (1991). 

[32] A. Aguirre and M. C. Johnson, Phys. Rev. D72, 103525 (2005), gr-qc/0508093. 

[33] J. J. Blanco-Pillado, M. Bucher, S. Ghassemi, and F. Glanois, Phys. Rev. D69, 103515 (2004), hep-th/0306151. 

[34] C. D. Carone and A. H. Guth, Phys. Rev. D42, 2446 (1990). 

[35] W. Israel, Nuovo Cim. B44S10, 1 (1966). 

[36] A. Vilenkin, Phys. Lett. B133, 177 (1983). 

[37] J. Ipser and P. Sikivie, Phys. Rev. D30, 712 (1984). 

[38] D. Langlois, K.-i. Maeda, and D. Wands, Phys. Rev. Lett. 88, 181301 (2002), gr-qc/0111013. 
[39] D. H. Lyth, Phys. Rev. Lett. 78, 1861 (1997), hep-ph/9606387. 

[40] M. W. Choptuik, Physical Review Letters 70, 9 (1993). 

[41] A. Aguirre, T. Banks, and M. Johnson, JHEP 08, 065 (2006), hep-th/0603107. 
[42] T. Banks and M. Johnson (2005), hep-th/0512141. 
[43] T. Banks (2002), hep-th/0211160. 

[44] X. Dong and D. Harlow, JCAP 1111, 044 (2011), 1109.0011. 

[45] G. Calabrese, L. Lehner, O. Reula, O. Sarbach, and M. Tiglio, Class.Quant.Grav. 21, 5735 (2004), gr-qc/0308007. 
[46] L. Lehner, D. Neilsen, O. Reula, and M. Tiglio, Class.Quant.Grav. 21, 5819 (2004), gr-qc/0406116. 
[47] L. F. Richardson, Phil. Trans. Roy Soc. A 210 (1911). 

[48] E. J. Copeland, M. Glciscr, and H.-R. MuUer, Phys.Rev. D52, 1920 (1995), hep-ph/9503217. 

[49] G. Fodor, P. Forgacs, and M. Mczei, Phys.Rev. D81, 064029 (2010), 0912.5351. 

[50] M. A. Amin, R. Easther, H. Finkel, R. Flauger, and M. P. Hertzberg (2011), 1106.3335. 

[51] M. A. Amin, R. Easther, and H. Finkel, JCAP 1012, 001 (2010), 1009.2505. 

[52] M. A. Amin (2010), 1006.3075. 

[53] M. A. Amin and D. Shirokoff, Phys.Rev. D81, 085045 (2010), 1002.3380. 

[54] I. Dymnikova, L. Koziel, M. Khlopov, and S. Rubin, Grav. Cosmol. 6, 311 (2000), hep-th/0010120. 

[55] M. Hindmarsh and P. Salmi, Phys.Rev. D77, 105025 (2008), 0712.0614. 

[56] M. Khlopov, R. KonopUch, S. Rubin, and A. Sakharov (1998), hep-ph/9807343. 



37 



[57] D. Baumann, A. Dymarsky, I. R. Klebanov, L. McAllister, and P. J. Steinhardt, Phys. Rev. Lett. 99, 141601 (2007), 
0705.3837. 

[58] A. Liiidc and A. Wcstphal, ,JCAP 0803, 005 (2008), 0712.1610. 

[59] L. Boubekeur and D. H. Lyth, JCAP 0507, 010 (2005), liep-pli/0502047. 

[60] M. W. Choptuik, Experiences with an adaptive mesh refinement algorithm in numerical relativity. (1989), pp. 206-221. 
[61] S. J. Parke, Phys. Lett. B121, 313 (1983). 
[62] K. Copsey (2011), 1108.2255. 
[63] K. Copsey (2011), 1109.4931. 



